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Abstract. We present a model for the Lyman-a ab- 
sorbers that treats all objects (from the low-density for- 
est clouds to the dense damped systems) in a unified 
description. This approach is consistent with an earlier 
model of galaxies (luminosity function, metallicity) but 
also with the known description of the density field in the 
small-scale non-linear regime. We consider two cosmolog- 
ical models: a critical universe = 1 with a CDM power- 
spectrum, and an open CDM universe with f2o = 0.3, 
A = 0. We reproduce the available data on column den- 
sity distribution as a function of redshift, the value of the 
main new parameter, the background ionizing UV flux, 
being consistent with the observed limits. This allows a 
quantitatively trustable analytical description of the opac- 
ity, mass, size, velocity dispersion and metallicity of these 
absorbers, over a range of column densities spanning 10 
orders of magnitude. Moreover, together with an earlier 
model of galaxy formation this draws a unified picture of 
the redshift evolution of structures in the universe, from 
underdense clouds to massive high density galaxies, from 
weak to very deep potential wells. 



In this paper, we seek a quantitative analytical de- 
scription of the number of these objects as a function of 
their column density, and of their internal properties. For 
long, the only tractable analytical approximation to de- 
scribe gravitational condensations was the Press-Schechter 
(1974) approach that relates the number of condensed ob- 
jects to the early, linear, density fluctuations. This ap- 
proximation, however, is not quite suited to our purpose 
since in principle it only describes objects with a density 
contrast that has just reached the virialization threshold 
A c ~ 200. Our aim, thus, is beyond the reach of this 
approximation. Indeed, we have undertaken this task to 
take benefit of the recent progress (Valageas and Schaef- 
fer 1997, VS I) that allows one to describe matter con- 
densations of any density contrast directly in terms of 
the (non-linear) density field at the epoch under consid- 
eration, the latter in turn being related in a known way 
(Balian & Schaeffer 1989, Bernardeau 1994, Bouchet et 
al.1991, Colombi et al.1997) to the initial conditions. This 
approach has already successfully reproduced the lumi- 
nosity function of galaxies (Valageas and Schaeffer 1998, 
VS II). 



Key words: large-scale structure of Universe - Galax- 
ies: evolution - intergalactic medium - quasars: absorption 
lines 



1. Introduction 

Lyman— a absorption lines along the line-of-sight of re- 
mote quasars are due to a vast class of different objects, 
ranging from small subgalactic, very underdense (down 
to 1% of the average density of the universe!) still ex- 
panding hydrogen clouds up to the halos of very large 
and very overdense (10 4 above the mean) galaxies. Some 
have reached virial equilibrium, others are UV heated and 
strongly coupled to their environment. 



2. Multiplicity functions 

As was done for the galaxy multiplicity function in VS II, 
we consider a cloud of mass M at redshift z to be char- 
acterized by a density threshold A(M, z) that is defined 
once the physical properties of the Lyman— a clouds are 
identified. It will turn out to be more convenient to attach 
to each such cloud a parameter x given by 

1 + A(M, z) 



x(M, z) = = 



t[R(M,z),z] 



(1) 



where 



d 3 ri d 3 



r-2 



V 2 



6(ri,r 2 ) with 



is the average of the two-body correlation function 
£2(1*1^2) over a spherical cell of radius R and provides 
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the measure of the (highly non-gaussian since we consider 
the actual density field) fluctuations in such a cell. The 
function A(M, z) defines the external boundary R(M, z) 
of the clouds, and is not necessarily given by the viri- 
alization constraint A = A c . It will be specified in the 
next section. Then we write for the multiplicity function 
of these clouds at a given rcdshift z (see VS I) , in physical 
coordinates: 



, , , . dM p , . , , 
V(M, z)j^ = Jjxh{x)dx 



(2) 



where p is the mean physical density of the universe at 
rcdshift z, while the mass fraction in halos of mass between 
M and M + dM is: 

fj,(M, z)j^- = xh(x)dx (3) 

For later convenience, it will be useful to define 

r 1 {x,z) d ^ = J JL_xh{x)dx (4) 

as the physical number density of these halos at rcdshift 
z, the mass M(x, z) being specified by the choice of the 
density contrast A(M, z). These multiplicity functions de- 
scribe the number density of objects of scale M — M+dM, 
defined by the constraint A(M, z). 

The scaling function h(x) is a universal function that 
depends only on the initial spectrum of fluctuations, and 
has to be taken from numerical simulations although its 
qualitative behaviour is well-known since it bears some 
general model-independent properties whose origin is the- 
oretical (Balian & Schaeffer 1989) and which have been 
well checked against observations as well as simulations 
(Valageas et al.1999; Colombi ct al.1997): 

x < 1 : h{x) oc x u - 2 , x > 1 : h(x) cx x^^ 1 e~ x / x ' 



with to ~ 0.5, ui s 



-3/2, 



10 to 20 and 



f 

Jo 



x h(x) dx = 1 



(5) 



Note that h(x) has to be measured only once, for a unique 
length scale and epoch provided it is in the highly non- 
linear regime. Then the scale-invariance of the many-body 
correlation functions which is the basis of this approach 
(see VS I) allows one to derive the multiplicity functions 
for any scale and time in the non-linear regime, as can be 
seen from (4). The correlation function £, that measures 
the non-linear fluctuations at scale R can be modelled in 
a way that accurately follows the numerical simulations 
(see VS I for more details). We first consider the case of 
a critical universe = 1 with a CDM power-spectrum 
(Davis et al.1985) normalized to erg = 0.5. As usual we 
define erg as the value of the amplitude of the density fluc- 
tuations at scale 8/i _1 Mpc given by the linear theory 



(a 2 ^ £). We use the scaling function h(x) obtained by 
Bouchet et al.(1991). We choose a baryonic density pa- 
rameter fib = 0.04 and Hq — 60 km/s/Mpc. In fact, a 
power-spectrum P(k) oc k n with n = — 2 and the same 
normalization gives very similar results. Then we study 
an open CDM universe, with erg = 0.77, fib = 0.03 and 
Ho = 60 km/s. These values are those we used previously 
to build a model for galaxy formation and evolution (see 
VS II), so that we obtain an overall consistent picture of 
the universe over a large range of object masses and scales. 

We can note that, as will be clear from our results, 
the properties of the Lyman-a clouds we shall obtain 
depend mainly on the physical model we build to de- 
scribe these objects (that is how one defines, or recognizes, 
such absorbers) and not much on the detailed mass func- 
tion of gravitational structures, provided the above model- 
independent properties are fulfilled. The latter is of course 
necessary, in order to count halos which lead to various 
absorption features and to make sure that there is no in- 
ternal inconsistency: one must not count the same mass 
several times while keeping track of all the mass in the 
universe (even if not all the mass produces Lyman-a ab- 
sorption lines). This is not an obvious task since one needs 
to simultaneously count different types of objects defined 
by several criteria. The formulation (4) allows one to per- 
form such counts in a consistent way. In our model, a 
given column density can usually be produced by many 
different clouds and the number of absorption lines of a 
given equivalent width will depend on an integral of the 
mass function (with a suitable weight) over a large range 
of parent halos. This will lessen the dependence on the 
exact slope of h(x). Hence our results are certainly very 
general and robust, provided the physical picture we de- 
velop in this article is correct. This seems confirmed by 
the good agreement we obtain with observations, and this 
would hold as well for any power-spectrum not too differ- 
ent from CDM or P(k) oc fc~ 2 . 



3. Properties of Lyman-a clouds 

3.1. Small low-density clouds: Lyman-a forest 

We assume that the IGM in small low-density halos is 
heated by the UV background radiation J{v) to a tem- 
perature T = 3 x 10 4 K. As a consequence, baryonic den- 
sity fluctuations are erased in low-density regions or halos 
with a virial temperature T < T over scales of order Rd 
with: 



Rd{z) 



1 



1 



fim p 



(6) 



where C s is the sound speed, jc s ~ 5/3, t H the age of the 
universe and m p the proton mass. For density contrasts 
of order unity Rd is also the usual Jeans length. Hence, 
we consider a first class of objects of typical scale Rd, 
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characterized by their mass M, and such that their asso- 
ciated gravitational potential, measured by the associated 
virial temperature T„, r , is lower than To (while their ac- 
tual gas temperature, due to UV heating, is To). This de- 
fines the average density ph and density contrast A(M, z) 
of these absorbers. By definition the baryonic density is 
nearly constant over their extent. These patches of matter 
can be regions which are still in expansion, or even un- 
derdense areas, and small virializcd halos. One can note 
that an alternative approach to this "smoothing" of the 
baryonic density field on scales smaller than the "Jeans 
length" could be to suppress the high fc-modes of the ini- 
tial density field, for instance Bi & Davidsen (1997) used 
PiGM(k) = Pdm (fc) [1 + (k/kj) 2 ]~ 2 . However, this is not 
possible in our present study since we intend to model si- 
multaneously, in a consistent way, all Lyman-a absorbers. 
Indeed, in our model, Lyman limit and damped Lyman 
systems (which we shall describe in the next sections) are 
high density halos with a radius and an impact parameter 
which is usually smaller than the length scale Rd but with 
a large virial temperature. Hence, small scale density fluc- 
tuations play a crucial role for these high column density 
objects and must be properly taken into account. 

Note that although we describe here the Lyman-a for- 
est lines in terms of distinct "objects" of size Rd, we would 
obtain the same results by directly considering the fluctu- 
ations of the density field. Indeed, the intersecting length 
of the line of sight with a filament is typically of order Rd, 
even if the overall length of this extended object could be 
much larger. This is due to the facts that i) density fluc- 
tuations are larger on smaller scales and ii) that such a 
filament would not appear as a straight line but rather as 
a random walk. In addition, we note that we actually take 
into account all the matter which is distributed over the 
line of sight. These points are discussed in more details 
in the Appendix. Assuming that photo-ionization equilib- 
rium is achieved, the neutral number density nni in the 
halo is: 



riHi 



«(r ) 

Gi J21 



(i-y)(i-y/2) 



f2o "*p 



(7) 



where a(T ) = 4.36 lO" 10 ^ -0 '^ s" 1 cm 3 is the re- 
combination rate, and G\ J21 is the ionization rate of 
neutral hydrogen (see Black 1981). As usual, J 2 i is the 
UV background radiation at the HI ionization threshold 
(912A) in units of 1CP 21 erg s _1 Hz -1 cm~ 2 sr -1 , and 
G 1 = 3.2 1CT 12 (Haardt & Madau 1996). Here n b /fl is 
the ratio of the baryonic density to the total density and 
Y = 0.26 is the helium mass fraction. Using the parameter 
x, we write 



tihi = ni(l + A) 2 = m x 2 £(R d ) 2 



where we define: 
a(T ) 



m(z) = 



Gi J 21 (z) 



(l-T)(l-F/2) 



Qb_p{z) 



(8) 



(9) 



As required by the Gunn-Peterson test, a given baryon 
mass fraction is quite inefficient in producing neutral ab- 
sorbing gas: tihi <C rib, for small ovcrdensitics. Since by 
definition the baryonic density is roughly constant over 
the whole region, we neglect the influence of the impact 
parameter of the line of sight which intersects the con- 
sidered patch of matter, and we assign to this region a 
constant neutral column density Nhi- 



Nhi = 2 riHi 2/3R d = 4/3 m x 2 £(i? d ) 2 R d 



(10) 



(the factor 2/3 comes simply from the average over the 
lines of sight of the depth of a spherical cloud). The num- 
ber of such regions which a line of sight intersects per 
redshift interval is: 



dn = irR 2 , c^-dz riix, z) — 
dz x 

Using (4) we obtain: 



d 2 n 



dlnNnidz 



3 1 dt 
-c- 



8R d dz£(R d ) 



xh(x) 



(11) 



(12) 



Thus, the slope of the column density distribution 
dn/dlnNni depends on the scaling function h(x), or more 
generally on the multiplicity function of mass condensa- 
tions. As we shall see below, these "objects" are low den- 
sity regions A ~ 1, so that i<1. Hence we are in the do- 
main where h(x) cx x u ~ 2 , which together with (10) leads 
to: 



dlnNnidz 



cx N 



(u>-l)/2 
HI 



(13) 



with oj ~ 0.3 for n ~ — 2 (Colombi ct al.1997) where 
n is the slope of the initial power-spectrum P(k) cx k n 
(generally speaking w is a function of n). For very low 
density regions, however, the multiplicity function is no 
longer given by the power-law tail of h(x) and shows a 

cutoff for (1 + A) ~ £ i * J ^ 1 -C 1, corresponding to 
very underdense objects surrounded by regions of even 
lower density. This implies a lower cutoff for the column 
density distribution at 



N 



Hlllower 



= 4/3 m R d W 2aV(1 ~ w) 



(14) 



Note that although some of the clouds described in this 
regime are underdense (A < 0), they correspond to well- 
defined objects. Indeed, in the non-linear regime (that 
is on small scales or at late times), most of the vol- 
ume of the universe is formed by very low-density areas 

— — cj/(1— u;) 

p ~ £ />«/) that are even less dense than these 

clouds. In fact, on these deeply non-linear scales, the av- 
erage density p loses the significance it has on large scales, 
in the sense that it does not define any longer a density 
boundary between two physically different classes of ob- 
jects. 
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On the high column density side, this description of 
the hydrogen clouds is valid until the virial temperature 
T reaches T . Thus it stops at: 



1 + A12 = 45 



or 



Xl2 



45 



(15) 



which corresponds to an upper column density cutoff: 



Nniupper = 2700 m R d 



(16) 



In the case of a critical universe, with a power-spectrum 
index n ~ — 2 which implies (Colombi et al.1996) lo ~ 0.3, 
we obtain for a constant UV background: 

R d (z) = 96(1 + z)- 3/2 kpc , A12 = 44 



a:i2 = 0.05(1 + zy-' , N HIlower = 2 10 b (l + z) b ' s cm 
iVff/upper = 1.5 10 13 (1 + z) 4 - 5 cm" 2 



(17) 



These regimes indeed correspond to low-density halos. The 
low column density cutoff Nuiiower is usually too low to 
be observed. The upper cutoff N HIupper increases strongly 
with rcdshift, because the HI density rim ~ ni(z) is 
proportional to the square of the baryonic density which 
varies as p b oc (1 + z) 3 . We note that although the den- 
sity contrast of these regions is not very large (A < 44), 
£ 12 ~ 1000 (l+z) -2 - 7 is significantly larger than unity: one 
is actually in the deeply non-linear regime where the den- 
sity field has been greatly distorted. Thus, regions with a 
small density contrast may in fact have undergone the in- 
fluence of strong non-linear effects and cannot be modeled 
by linear theory or methods relying on simple corrections 
to the latter. For instance, it is likely that some clouds 
with a (dark matter) density close to, or even lower than, 
the average density of the universe p have in fact already 
"collapsed" (or been deeply disturbed by shell-crossing) 
and shocked their baryonic content. This would also af- 
fect the gas temperature. 

3.2. Galactic halos 

Massive halos with a virial temperature above T do not 
see their baryonic density profile smoothed out via the 
heating produced by the UV background. Hence they con- 
stitute a second family of patches of matter, where we have 
to take into account the variation of the density with the 
impact parameter of the line of sight. Thus, we consider 
that virialized halos have a mean density profile p oc r~ 7 
with 7 = 1.8, which is consistent with the flat slope of 
the circular velocity observed in spiral galaxies, as well as 
with the observed galaxy correlation function. Indeed, for 
deep potential wells we expect 7 to correspond to the slope 
of the correlation function £. In fact, massive galactic ha- 
los must satisfy simultaneously the virialization constraint 
and a cooling condition (Silk 1977, Rees & Ostriker 1977). 
Similarly to the model we developed in VS II for galax- 
ies, small halos with a low circular velocity verify that 



A = A c because their cooling time is small so that their 
radius is given by the virialization constraint, while mas- 
sive halos have a long cooling time so that their boundary 
is the cooling radius which is of the order of R coo i = 120 
kpc (this simply means that we associate Lyman-a clouds 
with galaxies and not with clusters of galaxies) . This de- 
fines the functions A(M,z), or R(M,z), which we intro- 
duced in Sect. 2. 

We use for the gas temperature T gas of these halos the 
prescription: 



T gas = Min ( T , 2 10 6 if ) 



(18) 



This means that below 2 10 6 K the temperature of a large 
part of the gas which can cool still remains of the order of 
the value it reaches through shock-heating while for higher 
temperatures cooling is so efficient it does not get much 
higher. This also ensures that we recover the temperature 
range obtained in numerical simulations (e.g. Miralda- 
Escude et al.1996). Note that the threshold 2 10 6 K also 
corresponds to the virial temperature where shock-heating 
due to supernovae stops playing a dominant role (see VS 
II). However, removing this upper cutoff for T gas leads to 
nearly identical results which shows that our model is not 
very sensitive to its exact value. Indeed, even at z <~ 
halos with a high virial temperature T > 2 10 6 K are 
rather rare so that most absorption lines come from shal- 
lower potential wells. We also model crudely the collapse 
of baryons within the dark matter halos by assuming that 
the gas is distributed within the potential well over a ra- 
dius smaller by a factor A < 1 than the dark matter radius 
with the same density profile p oc r~ 7 : 



R^XR 



(1 + A) -» A" 3 (1 + A) 



(19) 



We use A = 0.4 for deep potential wells with T > 5 10 4 
K and a continuous transition to A = 1 at T = Tq since 
for shallow objects T < T there is no collapse (the gas 
is photo-heated to a temperature larger than the virial 
temperature). 

The neutral number density at the radius r within such 
a halo of external radius R and density contrast A is: 



riHi = «o(l + A) 



/ r \ - 2 t 



\R 



where we defined no(z) by: 

n (z) = p?fLL(l-Y)(l-Y/2) 



(20) 



(l _ 1) ^fe pi 2 ) 1 2 

\ 3/ £1q mp 



The column density Njji along the line of sight which 
intersects a cloud of radius R at impact parameter b is: 



N HI = 2n R(l+A) 2 ( - 



I-27 f^(R/b) 2 



du 



(1+M 2 )7 



(21) 
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In order to simplify the calculations, we consider two limits 
in the previous relation: b <C R and b — > R. In the case of 
a small impact parameter we write: 



Note that by doing so we disregard the mass which is 
outside of the considered halo. The number of halos along 
the line of sight is now: 



N HI ~ 2n i?(l + A) 2 ( — 



1-27 



<9 3 r 



( 22 ) \dlnN HI dlnxdz) 



\ 3 1 di , 

= --=:C— (1 



2< 



2^ 1 + A ^ 



where we noted 1^ = / °°(1 + u 2 )~ 7 du. The number of 
halos of this kind which a line of sight intersects at the 
impact parameter 6 per redshift interval is: 



7 ™ 7 77 dt , , x dx 
an = zirbab c—az rj(x, z) — 
dz x 



(23) 



Thus we obtain for the contribution due to small impact 
parameters, 



d 3 n 



dlnNnidlnxdz 



1 dt 
c- 



2> 



2(27-1) R dz 



(1 + A) 



(5-2 7 )/(2 7 -l) 



N 



HI 



(2// 



-2/(27-1) 



(24) 



x 2 h(x) 



Here the index 2 refers to the fact that these clouds form 
our second population of objects (the Lyman-a forest de- 
scribed in the previous section is our first population) 
while > corresponds to the small impact parameter regime 
(hence high Njji for a given cloud). As was noticed by 
Rees (1986), in such a model the density profile of virial- 
ized halos p oc r -7 governs the slope of the distribution of 
column densities: 



8 2 n 



dlnNnidz 



oc TV 



-2/(27-1) 



HI 



(25) 



2> 



through the dependence of the column density produced 
by a given halo on the impact parameter b. This is con- 
sistent with observations for 7 ~ 2, which is indeed the 
case (for 7 = 1.8 we have a slope of —0.77). However, this 
power-law is only valid over a limited range in Nhi (which 
translates into the dependence on Nhi of the boundaries 
of the domain of integration in a;). Indeed, for a given 
halo, this regime only applies for impact parameters much 
smaller than the halo radius (so that (22) is valid) but 
larger than the critical radius R n where self-shielding be- 
comes important and hydrogen is mainly neutral. More- 
over, as we explained in the previous section, small halos 
cannot be described in this manner as the UV background 
smooths their density profile which introduces the addi- 
tional length scale R<t- 

In a similar fashion, when the impact parameter is very 
close to the radius of the halo, we write: 



b = 



i-i^Va+A)-* 

2 \2n Q RJ V ' 



(27) 



N 



HI 



2n R 



x 2 h{x) 



Here the index < refers to the large impact parameter 
regime. Thus, we obtain a completely different power-law 
for the column density distribution dn/dlnN NHI oc N HI , 
which does not depend on the halo density profile and is 
only due to geometrical effects. In fact, its precise form is 
not important and it mainly plays the role of a cutoff in 
the column density distribution: this simply means that a 
given cloud mainly produces column densities larger than 
a characteristic value. For a fixed column density Nhi 
we choose the transition between both regimes (small and 
large impact parameter) as the point where the numbers of 
halos d 3 n/dlnN H idlnxdz given by both approximations 
are equal. In term of the variable x it corresponds to 



x 2 : 1+A = ( 2 7 -l)( 27 - 1 )/^)/ oc 1 /(47) ( ElL ) (28) 
with 

(1 + A) = A" 3 (1 + A c ) if x < (1 + A c )/l(R cool ) 
. R = A R coo i if x > (1 + A c )/l(R cool ) 

As we can see in (22) and (26) it means b ~ R as it should: 
in fact one could simply use the formulae for b <C R up to 
b = R and stop there. 

3.3. Neutral cores 

Because of self-shielding, the deep cores of massive halos 
are not ionized: all photons of the external UV background 
are absorbed by the outer shells of the halo. Using (20), 
we define the optical depth t at the distance I from the 
center of the halo (which is assumed to be spherical) by: 



1/2 



/ R / r \ -27 

<r pl n (l + A) 2 {-) dr 



1-27 



r ~ ^—^a pl n n R(l + A) (^-J if I < R (29) 

where a P i is the photo- ionization cross-section. Thus, for 
each halo we define a "neutral radius" R n , where r = 1, 
which determines the extension of the neutral core: 



R 



(26) R n =R 



(T pi n R{\ + A)' 
2 7 -l 



1/(27-1) 



(30) 
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Within R n all the hydrogen is neutral. Hence nni is pro- 
portional to the baryonic density (and no longer to its 
square), and we obtain: 



(2n). This corresponds to b = R n , and using (22) to 
Nhi = N H in with 

N H in = 2(2 7 - 1)700(7"/ ~ 2 10 18 cm- 2 (37) 



r \ — 1 



n HI = n 2 (l + A) (1) 



with 

n 2 (z) = (l- 7 /3)(l-r)^^. 



(32) 



Then, we proceed exactly as we did for the ionized part 
of the halos in the previous sections. We again divide this 
population in two sub-classes according to the value of the 
impact parameter with respect to R n . For small impact 
parameters we have: 



<9 3 r 



dlnNnidlnxdz 



2n> 



2( 7 -1)R dz 



!Al + A) (3-7)/(7-l) 



(JL 



HI 



\2n 2 RJ a 



-2/(7-1) 



x 2 h(x) 



(33) 



where we defined Joe = J °°(l + u 2 )~' y / 2 du. Here the n 
in the index 2n > refers to the fact that we consider the 
neutral cores within the population (2) halos. As was the 
case in the previous sections for the ionized shells, the col- 
umn density distribution we obtain is a power-law (within 
a limited range) with a slope determined by the halo den- 
sity profile: 



d 2 n 



dlnNnidz 



oc N 



-2/(7-1) 
HI 



2n> 



The coefficient 2 7 is changed to 7 because within deep 
neutral cores the HI density is proportional to the bary- 
onic density nni oc pt, while in the ionized shells it is pro- 
portional to its square nni oc p\ in the regime of photo- 
ionization equilibrium. This leads to a slope steeper than 
found previously: for 7 = 1.8 we get —2.5 instead of —0.77. 
When the impact parameter is very close to R n we obtain: 



d A n 



2R dz { + ' 



dlnNnidlnxdz , 2n< 

a pl n U \ ( IS hi \ 



27- 1 



\2n 2 RJ 



x 2 h(x) (35) 



The transition between both regimes (small and large im- 
pact parameter) corresponds to: 



x n : 1 + A 



apiiioR 



7-1 



(Nhi \ 
27- 1 ) \2n 2 R J 



27-1 



which occurs for b <~ R n - 

Finally, we have to take care of the transition between 
ionized envelopes: regime (2), and neutral cores: regime 



(31) which gives indeed r = <T pi NHi/2 ~ 1. In order to 
take into account, in a crude way, the non-sphericity of 
clouds, we use a smooth transition between ionized and 
neutral regimes. Thus, we multiply the contribution (2) 
by a factor exp[—(NHi/NHi n — 1)A] where e is a pa- 
rameter of order unity which describes the irregularity 
of clouds. We shall use e = 5 in the subsequent calcu- 
lations. For Nhi < Nnin we do not change the value we 
obtained previously for the contribution (2), since the pos- 
sible non-sphericity of the cloud can only ionize the gas, 
which otherwise would be neutral, by decreasing the col- 
umn density along a few directions, while a cloud which 
would be ionized according to the spherical model remains 
ionized. In a similar fashion, we correct the contribution 
(2n) by a multiplicative factor smaller than unity, equal 
to 1 — cxp[—(NHn/NHin — l)/ 6 ], where Nhh is the col- 
umn density which would be obtained in the regime (2), 
if the cloud was considered to be ionized. It is worth not- 
ing that even with such a large spreading factor e there is 
little smearing out of the distribution of the column den- 
sities. In particular it is not sufficient to fill up the dip at 
]og(N HI ) ~ 19. 

3.4- Final column density distribution 

After this description of the different regimes which are 
involved in the Lyman— a absorption lines, we can obtain 
the number of clouds of a given neutral column density 
which a line of sight intersects per redshift interval. We 
simply have to collect all the contributions we developed 
above, as a given absorption line may be produced by 
different kinds of objects. The integrations over x of the 
various distribution functions d 3 n/ 1 dlnNnidlnxdz are de- 
limited by the range of validity of the physical regimes 
they represent, as seen previously. We also add a lower 
cutoff for the impact parameter b > R Cl in order not to 
count the galaxy luminous cores (this is because a line-of- 
sight towards a remote quasar is never chosen to cross the 
luminous core of a galaxy). We shall take in the numerical 
calculations R c (z) = 1 (1 + z)~ 3 / 2 kpc. In addition, the 
distribution functions d 3 n/ dlnNnidlnxdz allow us to get 
the distribution of a given absorption line (a fixed neutral 
column density Nhi) with mass, radius, impact parame- 
ter, or any other cloud property. This is not the case for 
clouds of the first class (1), where each column density 
corresponds to a specific cloud. Finally, we need to spec- 
ify the value of the UV background, and its evolution with 
redshift. We treat J 2 \{z) as a free parameter (see Tab.l), 
which we choose so as to reproduce the column density 
distribution observed at the relevant redshifts while being 
consistent with the observational constraints (in addition 
J 2 i(z) increases with time until z <~ 2 and then drops at 



(34) 



(36) 
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low redshifts, as predicted by usual models of structure 
formation, see for instance Haardt & Madau 1996). In- 
deed, we shall see in the next paragraph that we obtain 
the right shape for the column density distribution, and 
the UV flux J21 only gives the normalization. In fact, as 
is well known, the latter only depends on the combina- 
tion QI/J21 (apart from an additional temperature de- 
pendence) for Lyman forest clouds and Lyman limit sys- 
tems, which is clearly seen on the definitions of n$(z) and 
ni(z), see (9). However, this is no longer true for damped 
systems, which consist of neutral cores embedded within 
massive halos, and only depend on f2b, see (32). 

4. Numerical results, J~2 = 1 

We first consider the cloud properties we obtain with our 
model in a critical universe £1 = 1 with a CDM power- 
spectrum (Peebles 1982, Davis et al.1985) normalized to 
cr 8 = 0.5. 

4-1. Redshift evolution of the column density distribution 

Fig.l shows the contributions (1), (2) and (2n) to the col- 
umn density distribution, at z = 0, 2.5 and 5. We can see 
very clearly on this figure the 3 classes of objects that we 
defined in Sect. 3. 
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Fig. 1. The contributions (1) corresponding to the Lyman-a 
forest (top left), (2) corresponding to the Lyman limit systems 
(central part) and (2n) corresponding to the damped Lyman 
systems (bottom right), to the column density distribution, at 
z = (solid line), z = 2.5 (dashed line) and z — 5 (dot-dashed 
line) . 

Low column density lines (Nhi < 10 15 cm -2 at z = 2.5 
and N HI < 10 16 cm" 2 at z = 5, on the left of the picture) 
correspond to low density areas and shallow potential wells 
(contribution (1)) with a roughly uniform density on scales 



smaller than the "damping" length Rd- The importance of 
this contribution decreases at small z and shifts to low col- 
umn densities because as time goes on most of the matter 
becomes embedded in deep potential wells, which are de- 
scribed by other regimes, while the average density of the 
universe declines. This translates into the redshift depen- 
dence of the upper column density NHiupper(z), see (17). 
These objects exist down to quite small column densities 
(N HI > 10 10 cm" 2 at z = 2.5 and N HI > 10 12 cur 2 at 
z = 5), that is cloud dark matter masses as low as 1O 1O M0 
at z = and 10 6 M© at z = 2.5 and z = 5. The associated 
absorption lines correspond to the Lyman— a forest. 

Larger column densities, up to N HIn ~ 2 x 10 19 cm -2 , 
are described by the regime (2). They are produced by 
virialized halos, which the line of sight intersects with an 
impact parameter much smaller than their radius but large 
enough so that the hydrogen is highly ionized. The low col- 
umn density cutoff associated with this regime varies with 
redshift, as it corresponds to the upper limit of the regime 
(1) we described above, while the high column density cut- 
off is simply given by the ionization condition and is of the 
order of Njiin, which is the column density at which the 
optical thickness is unity and whence is independent of 
z. Note that the regime (2<) corresponding to the rising 
part on the left of contribution (2), see eq.(27), plays no 
role since in this domain most observed lines are produced 
by forest clouds (regime (1)) described above. Indeed, the 
transition (2>) - (2<) coincides with the transition (1) 
- (2) as it should: we switch (almost) continuously from 
one class of objects to the other one. The absorption lines 
described by this regime (2) correspond to Lyman limit 
systems. 

Finally, very high column densities (on the right of 
the figure), larger than Nijin, correspond to lines of sight 
which intersect neutral cores. There is a gap in the distri- 
bution of column densities because the fraction of neutral 
hydrogen switches suddenly from ~ 1CP 3 in ionized shells 
to 1 in neutral cores (this is linked to the sharp change of 
the factor exp(— r) which behaves as an exponential of the 
column density). As we explained in Sect. 3, we can see on 
the figure that the slope we obtain (for Nhi ~ 10 22 cm~ 2 ) 
is steeper than the one we found for regime (2) . These ab- 
sorption lines correspond to damped Lyman systems. 

We show in Fig. 2 the total column density distri- 
bution at z = 3 which results from the sum over all 
contributions. We can see that our model agrees rea- 
sonnably well with the observations, from Nhi ~ 10 13 
cm~ 2 up to Nhi ~ 10 22 cm -2 . In particular, we recover 
the under-abundance of lines at intermediate column den- 
sities Nhi ~ 10 15 — 10 16 cm -2 . This feature is also quite 
clear on Fig.l. It corresponds to the transition between 
populations (1) and (2). The gas embedded within galac- 
tic halos can cool and collapse (for deep potential wells 
T > T ) which increases the column density along the 
line-of-sight (the neutral number density scales as the 
square of the baryonic density). As a consequence, clouds 
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Fig. 2. The total column density distribution (summed over 
all contributions: Lyman-a forest, Lyman-limit and damped 
Lyman-a systems) at z = 3. The data points are from Hu et 
al.(1995) (crosses) and from Petitjean et al.(1993) (rectangles). 



to the gap in Fig.l. This is also seen in numerical simula- 
tions (Katz et al.1996). We can see that the slope of the 
distribution function we obtain is again consistent with 
observations. 



-2 



CDM , Q= 1 



N,„>1.6 10 17 



which would have produced column densities somewhat 
larger than Nniupper{z) actually lead to deeper absorp- 
tion features. This produces an under-abundance of ob- 
jects around Nhi ~ 10 15 cm" 2 at z = 2.5. 




Fig. 4. The evolution with redshift of the total column density 
cumulative distribution for Nhi > 10 14 cm" 2 (upper curve), 
Nhi > 1.6 10 17 cm" 2 (intermediary curve) and Nhi > 2 10 20 
cm" 2 (lower curve). The data points are from Bechtold (1994) 
(stars), Bahcall et al.(1996) (triangles), Giallongo et al.(1996) 
(filled squares), Kim et al.(1997) (open circles) and Lu et 
al.(1996) (filled circle) for Nhi > 10 14 cm" 2 ; from Stor- 
rie-Lombardi et al.(1994) (crosses) for Nhi > 1.6 10 17 
cm -2 ; and from Wolfe et al.(1995) (filled squares) and Stor- 
rie-Lombardi et al.(1995) (cross) for Nhi > 2 10 20 cm" 2 . 



Table 1. Redshift evolution of the background UV flux J2i(z) 
used in this article for both cosmologies (critical and open uni- 
verse) . 

7=0 1 2 3 4 5~ 

Q=l Jn = 0.05 0.4 0.7 0.4 0.2 0.1 
fio =0.3 0.05 0.5 0.8 0.4 0.2 0.1 



Fig. 3. The total column density cumulative distribution at 
z = (solid line), z = 2.5 (dotted line) and z = 5 (dot-dashed 
line). The data are from Bechtold (1994) (Forest at z = 2.5), 
Lanzetta et al.(1995) (Lyman limit systems at z = 2.6) and 
Wolfe et al.(1995) (damped Lyman systems at z = 2.25). 

Fig. 3 shows the total (summed over all contributions) 
column density cumulative distribution at z = 0, 2.5 and 
5. The flat part for 10 19 < Nhi < 10 20 cm" 2 corresponds 



The evolution with redshift of the total column density 
cumulative distribution is shown in Fig. 4 for Nhi > 10 14 
cm" 2 , N H i > 1.6 10 17 cm" 2 and N HI > 2 10 20 cm" 2 . 
We can see in the figure that we can fit the data simul- 
taneously for the 3 types of absorption lines which are 
representative of 3 different classes of objects: Lyman— a 
forest clouds, Lyman limit systems and damped Lyman 
systems, which are indeed described in our model by 3 dif- 
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ferent regimes. This builds our confidence in the validity 
of this model. Note that we could increase the number of 
damped systems by using a slightly larger o§ which does 
not significantly modify our results for limit and forest 
objects, as shown below in Sect. 4. 2. With = 0.04, the 
values of the UV flux we use are shown in Tab.l. They 
are consistent with observations: Giallongo et al.(1996) 
find J 2 i = 0.5 ± 0.1 for 1.7 < z < 4.1, while Cooke et 
al.(1997) get J 2 i = lto.l for 2 < z < 4.5, with no evi- 
dence for any redshift evolution within these intervals. At 
low redshifts z < 1 the UV flux shows a sharp drop (at 
z = Vogel et al.1995 find J 21 < 0.07 while Donahue et 
al.1995 get J 2 i < 0.033). Moreover, the overall redshift 
dependance of J 2 \{z) shown in Tab.l is similar to the pre- 
diction of the usual models of galaxy formation from hi- 
erarchical scenarios where the radiation comes from stars 
and quasars (e.g. Haardt & Madau 1996). We can see in 
Fig. 4 that we recover the observed break in the redshift 
evolution of the number density of lines at z <~ 1.5 (e.g. 
Jannuzi 1998). Indeed, at lower z the number of lines with 
Nhi > 10 14 cm~ 2 remains constant while it slightly in- 
creases for Lyman-limit and damped systems. This break 
in the continuous decline with time of the number of forest 
lines due to structure formation which builds increasingly 
deep potential wells is produced by the sudden drop of 
the UV flux at low z. At low redshifts z~0we may over- 
estimate the number of lines with Nhi > 1.6 10 17 cm~ 2 
and Nhi > 2 10 20 cm" 2 which come from galactic halos 
since we did not take into account star-formation which 
consumes (and may eject) some of the gas. 

We can note that we manage to get a satisfactory 
agreement with observations for the column density distri- 
bution of Lyman-a clouds, while using flf, = 0.04, Hq = 
60 km/s, cr 8 = 0.5 and J 21 (z = 2.5) = 0.7. On the 
other hand, numerical simulations need fi& > 0.05 for 
J 2 i(z = 2.5) > 0.2, cr 8 = 0.79, fl = 0.4, A = 0.6, and 
H = 65 km/s (Miralda-Escude et al.1996) or Q b > 0.05 
for J 2 i(z = 2.5) > 0.1, er 8 = 0.7, = 1 and H = 50 km/s 
(Katz et al.1996). These latter results would mean that it 
is difficult to satisfy the nucleosynthesis bounds with the 
observational estimates of J 2 i. The fact that we do not 
encounter such a serious problem here is rather encourag- 
ing and suggests, as we shall see below, that due to the 
very extended range of clouds (from weak potential wells 
to very deep halos) which contribute to a given Nhi it is 
difficult to take into account properly all contributions (es- 
pecially from the shallower clouds) in a simulation, which 
should thus usually underestimate the column density dis- 
tribution function. 

4-2. Influence of various parameters 

As we explained in Sect. 3, for all regimes the slope of the 
column density distribution only depends on the shape 
of the initial power-spectrum through ui and 7 (we take 
7 = 1.8 but it could slightly depend on the slope n of P{k)) 



as we can see in (13), (25) and (34). This constrains n to 
be close to —2 hence both a CDM-like power-spectrum 
(since the local slope on the scales of interest is indeed 
close to —2) and a power-law power-spectrum with n ~ 
—2 give satisfactory results. However, the normalization of 
the column density distribution depends on the UV flux 
J 2 i, as we noticed above, and on cosmological parameters. 
Indeed, we can see from Sect. 3 that we have: 

Kma^dzj^K-^-) URd)Rd] h 

(dlnN H idz) 2> ^ ( m ) 

K\dlnN H idzJ 2n> \ n 2 ) 
with 

rp-0.75(Q h 2\2 

rnoc^ ) b ' and n 2 ozn b h 2 (38) 

where we neglected the influence of the boundaries of the 
integrations over x. Fig. 5 shows the influence of various 
parameters on the column density distribution. If we in- 
crease the normalization o- 8 of the power-spectrum hence 
£ (dotted line in the figure), the number of Lyman-a forest 
lines decreases: clustering is more important (one is deeper 
into the non-linear regime) so that a larger fraction of the 
mass of the universe is within high density virializcd halos 
which produce Lyman limit and damped systems. Simul- 
taneously these two latter contributions increase, but this 
is not shown in (38) because we neglected the variation 
of the boundary in x. This increase is most important for 
highest column densities (massive damped systems) since 
these lines come from the rare high density halos which are 
in the exponential cutoff of the mass functions and whose 
abundance is very sensitive to the normalization of the 
power-spectrum (it enters as an exponential factor). On 
the other hand the number of Lyman-limit lines is nearly 
invariant because each of them is drawn from a large pop- 
ulation of possible parent halos: increasing o§ means there 
are fewer low mass halos but more large mass halos and 
both effects nearly cancel each other. If the baryonic den- 
sity parameter Qb gets higher (dot-dashed line) all con- 
tributions increase since there is more hydrogen available 
but the magnitude of this effect is not the same for all 
lines. Finally, if we increase J 2 i the number of lines will 
decrease in all regimes since the neutral fraction of hydro- 
gen is lower (dashed line). Once again, this is not seen in 
(38) for damped systems because it appears through the 
change of the boundaries over x: the deep cores of virial- 
ized halos are not influenced by a small increase of J 2i , 
which only destroys the outer neutral shells. Indeed we can 
check in Fig. 5 that the number of very large column den- 
sity lines does not change. We note that a change of Qb or 



10 



Valageas, Schaeffer & Silk: Lyman-a Clouds 



J21 simply leads to a horizontal translation of the column 
density distribution for forest and Lyman limit absorbers. 
Indeed, the number and properties (for the dark matter) of 
these halos do not change, but a given region will produce 
a column density which evolves with the baryonic fraction 
(fib) and the neutral hydrogen fraction (J2i)- This is no 
longer true for damped systems where a change of fib or 
J21 usually introduces new halos or removes the smallest 
ones. Moreover, we can see that the location of the gap be- 
tween Lyman limit and damped systems is constant since 
it only depends on the photo-ionization cross-section, see 
(37). 



which were already fully constrained by other sets of ob- 
servations. This implies for instance that the abundance of 
the largest damped lines is given by these previous mod- 
els. Indeed, as we explained above, these lines correspond 
simply to deep neutral shells of galactic halos, while the 
forest absorbers are new low-density objects (predicted by 
the same description of the non-linear density field, but in 
a different density regime) 

4-3. Mass and circular velocity associated to different 
clouds 
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Fig. 5. The column density distribution as in Fig.l at z = 2.5 
(solid line), with J21 = 3 instead of 0.7 (dashed line), with 
Qb = 0.1 instead of 0.04 (dot-dashed line), or with as = 1.3 
instead of 0.5 (dotted line), at the same redshift. 



Thus, we see that the normalization of the column den- 
sity distribution does not depend only on the usual param- 
eter ni, and that a change of a given variable (fib, J21, To 
or er 8 ) usually influences the three classes of objects de- 
scribed in Sect. 3 in a different way. The dependence on a s 
is very small for all values of interest except for the largest 
damped systems. This means that a power-spectrum nor- 
malized to the COBE data would also lead to a reason- 
able agreement with observations of Lyman-a clouds. The 
influence of fib and J21 is stronger, but it is degenerate 
for forest and Lyman limit lines through the combina- 
tion n\. Thus, for a given normalization of the power- 
spectrum, one could first derive fib from observations of 
large damped systems, and then obtain J 2 i from Lyman 
limit or forest absorbers, in order to match the data. How- 
ever, we must emphasize that we have only one important 
"free" parameter in our model: J 2 i (which must also be 
consistent with observations), since P(k) and fib are taken 
from consistent models of galaxies (VS II) and clusters 



Fig. 6. Evolution with redshift of the halo rotational velocity 
distribution function, for Nhi = 10 17 cm" 2 , at z — (solid 
line), z = 2.5 (dashed line) and z = 5 (dot-dashed line). 



Our model yields the evolution with redshift of the 
halo rotational velocity (or mass, radius,..) distribution 
function, for a fixed column density. We define the halo 
circular velocity v c as: 



GM 
R 



Gp(l + A)R 2 



(39) 



>l7 cm 2 at z 



and Fig. 6 presents the case Nhi — 10 
0, 2.5 and 5. The average halo rotational velocity (or mass) 
gets larger as time goes on, since gravitational clustering 
builds increasingly deep and massive potential wells. The 
sharp high velocity cutoff is not due to the multiplicity 
function of virialized halos but to the fact that very large 
and massive clouds cannot produce column densities as 
low as Nhi = 10 17 cm~ 2 . In other words, looking at a 
specific Nhi selects a finite range of parent halos (the 
contribution of larger clouds is not exactly zero because 
if the line of sight intersects such a halo very close to its 
external radius it can still produce a small column density, 
due to the small intersection length, but this occurs with 
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a negligible probability as can be seen in the figure). The 
location of this cutoff does not evolve much from z = 
to z = 2.5 because in this range clouds are defined by 
the cooling relation R — R CO oi- Indeed, for a given halo 
the characteristic "minimum" column density on a line of 
sight is given by: 



N 



HImin 



m(l + A) 2 i? 



RVQP 



n\R 



(40) 



similarly to the calculation described previously for regime 
(1) clouds, see (10). Hence, a given column density selects 
a "maximum" rotational velocity: 



J2i(z)N H iR; 



3 

cool 



1/4 



(41) 



which evolves slowly with z through JaiO^) • In fact, 
at low z this upper velocity cutoff increases with z be- 
cause of the sharp decline of 3i\{z) with time. Indeed, at 
z ~ 2 in order to observe a given column density (here 
Njji = f0 17 cm -2 ) one must look through a deeper po- 
tential well (larger mass and density) defined by a higher 
T and v c than at z ~ because J21 is much larger. On the 
other hand, at z = 5 the clouds corresponding to v cmax 
for N H i -- 
dition A 

Vcmax C< 



10 cm are defined by the virialization con- 
Ac, so that we get: 



which decreases strongly at high redshift (note that we 
neglected the influence of the collapse factor A). This be- 
haviour can be understood very simply within the frame- 
work of the physical picture we developed in Sect. 3 for 
Lyman-limit clouds, based on the association of these ab- 
sorbers with galaxies. Indeed, at high redshift as time goes 
on the hierarchical clustering process builds increasingly 
deeper and more massive but lower density potential wells 
(the density of these halos scales as (1 -I- A c )p), so that 
the rotational velocity attached to a given column den- 
sity gets larger with time, as in (42). However, after some 
time (typically z < 1 — 2) galaxies are no longer defined by 
the sole virialization constraint: massive galaxies are de- 
termined by a cooling constraint (VS II, Silk 1977, Rees & 
Ostriker 1977) which introduces a fixed length scale R coo i- 
As a consequence, the highest velocity associated with a 
given Njji does not evolve any more (except with the pos- 
sible changes of the background UV flux J 2 i) which leads 
to the regime described by (41). In other words, at small 
redshift we identify Lyman-limit clouds with galaxies, and 
not with clusters of galaxies. This leads to a qualitative 
change in the behaviour of some physical properties. Thus, 
the "simple" hierarchical clustering scenario (coupled with 
a proper physical picture) can produce non-trivial histo- 
ries, that is with qualitative discontinuities and amenable 
to observational checks. The lower cutoff is due to the 



fact that low velocity halos are described by regime (1): 
their density is roughly uniform over scales of the order of 
the "damping" length Rg and they cannot produce large 
column density absorption lines. Thus, it corresponds to 
the virial temperature To and shows (nearly) no evolu- 
tion with redshift. The small discontinuity in the curve at 
z = around v c — 60 km/s corresponds to the transition 
from halos defined by A = A c (virialization constraint) to 
R = Rcooi (cooling constraint). The latter discontinuity 
thus has no physical meaning and is simply due to our 
interpolation procedure between these different regimes. 
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Fig. 7. a) Evolution with redshift of the average rotational 
velocity of the halos associated to a given column density at 
2 = (solid line), 2 = 2.5 (dotted line) and 2 = 5 (dot-dashed 
line), b) Evolution with redshift of the average mass of the 
halos associated to a given column density, for the same cases. 



We note in Fig. 6 that the velocity distribution is quite 
extended, and that for z < 2.5 it is not far from uni- 
form from 20 km/s up to 400 km/s. This is due to the 
fact that a given cloud may produce many different ab- 
sorption features according to the value of the impact pa- 
rameter b. Conversely, a given column density Nhi may 
originate with a similar probability from many different 
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Fig. 8. a) Evolution with redshift of the mean halo rotational 
velocity, for clouds defined by Nhi = 10 14 cm -2 (solid line), 
Nhi = 10 17 cm" 2 (dashed line) and Nhi = 1.6 10 21 cm" 2 
(dot-dashed line), b) Evolution with redshift of the mean halo 
mass, for the same cases. 



clouds. Even at z = 0, small velocity clouds (v c < 100 
km/s) make up a sizeable fraction (~ 40%) of the total 
Lyman-a lines with Nhi = 10 cm" 2 . This could explain 
the fact that simulations with insufficient resolution pro- 
duce fewer Lyman-limit systems than is required by the 
observations (Katz et al.1996) since it appears that one 
should resolve clouds down to v c ~ 20 km/s. 

Alternatively, we can consider the evolution with red- 
shift of the mean halo rotational velocity or dark mat- 
ter mass associated with a given column density. This is 
displayed in Fig. 7 for z = 0,2.5 and 5. The feature at 
low velocity and mass corresponds to the transition be- 
tween regimes (1) and (2), while the sudden change at 
Nhi ~ 10 20 cm" 2 is due to the transition from lines of 
sight intersecting ionized shells to those probing the neu- 
tral cores of deep halos. For any column density the av- 
erage velocity and mass grow as time increases because 
the hierarchical clustering process builds deeper and more 
massive halos. At small Nhi in regime (1), there is a 



unique correspondance between column density and ve- 
locity, or mass, see (10), which gives: 

Nhi oc (1 + A) 2 oc v 2 oc M 2 (43) 

since the radius R = Rd is constant. For larger Nhi there 
is no longer such a unique relation, as different clouds can 
produce the same column density. The mean velocity or 
mass first increases with Nhi and then reaches a plateau 
where it is nearly independent of Nhi- The rising part 
of the curve corresponds to the fact that larger column 
densities can be produced by deeper and more massive 
halos (as we described above for the high cutoff of the ve- 
locity distribution at fixed Nhi)- Hence, as one looks for 
larger equivalent widths one adds to the population of par- 
ent halos new more massive and deeper clouds (while the 
minimum mass of the possible halos does not change since 
it is given by the fixed transition to the regime (1), one 
simply needs to draw lines of sight which pass closer to the 
center of this small potential well). As a consequence the 
average velocity (mass) increases with Nhi , until the high 
velocity (mass) cutoff gets larger than the typical velocity 
(mass) of objects which have collapsed at the considered 
redshift (this is where the cutoff of the mass functions 
comes in). Then, adding a new population of more mas- 
sive halos when looking at a larger Nhi only produces a 
negligible change in the mean velocity or mass because the 
number of these new clouds is very small. 

Finally, there is a transition at Nhi ~ 10 20 cm" 2 
where the average velocity and mass suddenly decrease 
to grow again later with Nhi- Indeed, for column densi- 
ties slightly larger than 10 20 cm -2 we start to probe the 
neutral cores of virialized halos. The sharp change in the 
fraction of neutral hydrogen reflects that column densi- 
ties Nhi ~ 10 20 — 10 21 cm" 2 can be produced by small 
clouds which were characteristic of the lowest Lyman-limit 
lines: the maximum velocity or mass allowed for a cloud 
to produce such a line decreases suddenly to the value 
corresponding to the transition between regimes (1) and 
(2), hence we recover the same mean velocity and mass, 
as can be seen in the figure. Then, exactly as occured 
for Lyman-limit systems, the average velocity (mass) first 
increases with Nhi and finally reaches a plateau, due 
to the cutoff of the mass functions. This final velocity 
(mass) is larger than the one reached in the regime (2) 
for Lyman-limit lines, because the velocity (mass) dis- 
tribution is more heavily weighted towards large velocity 
(mass). This is due to two factors. First, the neutral hy- 
drogen density nni is now proportional to the baryonic 
density and not to its square, see (20) and (31), which 
means that for a given cloud the change with the impact 
parameter b of the column density is slower but also that 
for slightly larger clouds the impact parameter must in- 
crease faster in order to produce the same Nhi, which 
leads to a cross-section factor b 2 more heavily weighted 
towards large clouds. Second, the neutral radius R n grows 
faster than R, see (30), which again favors massive halos. 
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Of course, this translates directly into the R and A de- 
pendence of d 3 n / dlnN h idlnxd z in both regimes, see (24) 
and (33). 

Fig. 8 shows directly the redshift evolution of the aver- 
age velocity or mass for 3 types of clouds, representative of 
the 3 classes of objects we described previously. The main 
feature is that for all clouds the mean velocity and mass 
decrease at higher redshift, since only weaker and smaller 
potential wells were formed at earlier times. This is not 
always true for the high column density Njji = 1-6 10 21 
cm~ 2 because of the transition between ionized shells and 
neutral cores which introduces additional effects, as de- 
scribed above. This also explains why the average mass or 
velocity associated with this damped system is smaller at 
low z than those characteristic of a lower column density 
Lyman-limit system. 

4-4- Radius and impact parameter associated to various 
column densities 
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Fig. 9. Evolution with redshift of the halo radius distribution 
function, for Nhi = 10 17 cm -2 , at z = (solid line), z — 2.5 
(dotted line) and z — 5 (long dashed line) . The filled rectangle 
and the small ellipse represent the contribution of the clouds 
defined by the constraint R = R coo i at z = and z — 2.5, see 
text. 

We present in Fig. 9 the redshift evolution of the halo 
radius distribution function for a fixed column density 
Njji = 10 17 cm -2 (here R is the radius of the dark mat- 
ter halo, i.e. without collapse). Note that massive halos are 
defined by a constant radius R = R coo i- In order to take 
into account this population in the picture, we also plot- 
ted in the figure the point (R coo i, d 2 n/ dlnNjjidz) due to 
this contribution shown as a filled rectangle for z = and 
as a small ellipse for z = 2.5. Thus, while the total con- 
tribution of halos defined by the virialization constraint 
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Fig. 10. a) Evolution with redshift of the average radius and 
impact parameter of the halos associated to a given column 
density at z — (solid line), z — 2.5 (dotted line) and z = 5 
(dot-dashed line), b) Evolution with redshift of the average 
neutral fraction of the halos associated to a given column den- 
sity, for the same cases. 

A = A c is given by the integral over InNjji of the curves 
shown in the figure, the total contribution of halos de- 
fined by the cooling constraint R = R CO oi is given by the 
ordinate of the point {R coo i,d 2 n/dlnNjjjdz), multiplied 
by a Dirac S(R — R CO oi) function, that in our more real- 
istic model (VS II) is somewhat smeared out. The latter 
does not appear for z = 5 because it is negligible. This 
allows us to get at once the proportion in number of both 
classes of objects. Thus, one can see that at z = the 
"Rcooi population" dominates so that most clouds asso- 
ciated with Nhi — 10 17 cm -2 have a dark-matter halo 
radius R ~ R CO oi = 120 kpc. At higher redshift the char- 
acteristic radius decreases but one can see that at z = 2.5 
the radius distribution is still quite extended, ranging from 
4 kpc up to 120 kpc. The low radius cutoff corresponds to 
the transition with the regime (1) or to the influence of 
the cutoff radius R cu t (indeed to observe a high column 
density through a small halo one should look close to the 
center where the gas density is high, however this is not 
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Fig. 11. a) Evolution with redshift of the mean radius, for 
clouds defined by N H i = 10 14 cm -2 (solid line), N H i = 10 17 
cm -2 (dashed line) and Nhi = 1.6 10 21 cm" 2 (long dashed 
line), b) Evolution with redshift of the mean impact parameter, 
for Nhi = 10 17 cm -2 (dashed line) and Nhi = 1.6 10 21 cm" 2 
(dot-dashed line). 

always allowed since we remove lines of sight which would 
cross the galactic luminous core). It grows with time in 
parallel with the "damping" length Rd and the lower cut- 
off R c . Note that there is also in our present model an 
accumulation of clouds with radius Rd (but various col- 
umn densities) corresponding to population (1). 

Fig. 10 shows the average radius and impact parame- 
ter of the halos associated with a given column density at 
z = 0,2.5 and 5. The characteristic radius and impact pa- 
rameter decrease at higher z since virialized objects were 
smaller in the past. The panel a) shows the same features 
for the average radius of Lyman- limit and damped systems 
as those displayed in Fig. 7 for the circular velocity or mass, 
that is at first increasing with Nhi followed by a plateau, 
for the same reasons as those described previously. For 
Lyman forest clouds the radius R = Rd is constant by 
definition. As shown in panel b), the impact parameter 
decreases for higher column densities, both for Lyman- 
limit and damped systems, since one has to probe deeper 



towards the center of the clouds to reach sufficiently high 
neutral hydrogen densities in order to obtain large Nhi 
(the temporary increase within the range corresponding 
to Lyman-limit systems is due to the collapse of baryons 
which leads to larger b for a given column density because 
the gas density is larger). The feature at Nhi ~ 10 20 — 10 21 
cm" 2 is due to the transition between ionized shells and 
neutral cores, as usual. The impact parameter is not dis- 
played for forest clouds (regime (1)) since it plays no role 
for this population and was not specifically defined (the 
radius is the sole relevant scale). We note that Bechtold 
et al.(1994) derive from observations at z = 1.8 a radius 
57 < R < 400 kpc, with a median value of R = 130 kpc, 
for clouds with N H i ~ 10 14 cm" 2 . Dinshaw et al.(1994) 
obtain similar values. The average radius we obtain in 
Fig. 10 is smaller (~ 21 kpc) but as shown in Fig. 9 the 
radius distribution function is quite extended so that one 
indeed expects to observe clouds with a radius up to ~ 100 
kpc. Moreover, as discussed in Appendix, our modelling is 
compatible with quite elongated filamentary clouds. The 
radius we refer to should be considered as the average in- 
tercept along the line-of-sight. However, if these very un- 
derdense clouds are filamentary, the distance at which two 
separate lines can hit the same cloud is much larger. So, 
we rather consider this offset as a possible manifestation 
of the non-sphericity of the absorption features. 

Fig. 11 shows the redshift evolution of the average ra- 
dius and impact parameter for 3 column densities. One can 
see clearly the decrease at high z of both length scales. 
The decline of b is slower because at larger redshift the 
UV flux is lower while characteristic densities are higher, 
which tends to increase the impact parameter. We can also 
note that very different Nhi have very close mean radii 
and impact parameters, while the radius and impact pa- 
rameter distributions for a fixed column density are very 
extended, as was shown in Fig. 9. Of course, this is due 
to the fact that a given Nhi can be produced by a large 
range of clouds, so that very different column densities are 
in fact drawn from the same population of halos. 

4-. 5. Metallicity 

Fig. 12 shows the average metallicity (that is the abun- 
dance of oxygen or any other element mainly produced by 
SN II) of the halos associated with a given column den- 
sity at z = 0,2.5 and 5. We obtain this metallicity from 
our model for the evolution of galaxies described in VS II. 
The latter is entirely consistent with the present study of 
Lyman— a absorption lines, as it uses the same descrip- 
tion for the multiplicity functions of mass condensations. 
It also involves the same model of star-formation which en- 
abled us (VS II) to get the luminosities of galaxies, as well 
as their metallicity, with no new parameter. In fact, within 
this framework we define three metallicities, correspond- 
ing to stars (Z s ), star-forming gas concentrated within the 
inner parts of galaxies (Z c ), and diffuse gas spread over 
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Fig. 12. The average metallicity Zh of the halos associated to 
a given column density at z — (solid line), z — 2.5 (dashed 
line) and z = 5 (dot-dashed line). 



the halo (Zh) (corresponding to the two-component model 
of VS II). Lyman- limit lines, which arise when the line of 
sight intersects the outer ionized shells of a galactic halos, 
are characterized by Z^, while for damped systems, corre- 
sponding to neutral cores, we should observe a metallicity 
in the range Zh to Z c . We do not assign a metallicity to low 
column density forest lines, corresponding to regime (1), 
as they are associated with low density regions which may 
have not virialized yet. Observations by Lu et al.(1998) 
also find that there is a sharp drop in the metallicity of the 
gas (although for C) from [C/H] ~ -2.5 for N m > 10 14 ' 5 
cm~ 2 downto [C/H] < -3.5 for N HI < 10 14 cm" 2 , at 
redshifts 2.2 < z < 3.6. Note that we predict a significant 
increase with z of the column density corresponding with 
this transition. One may expect that diffusion processes 
(through winds or ejection of gas during violent mergers) 
would slightly lower the column density associated with 
this drop and would produce a non-zero metallicity for 
objects described by regime (1), which may also be en- 
riched by a generation of population III stars. Naturally, 
the metallicity decreases at high redshift when star forma- 
tion has not had enough time to synthesize many metals. 
At a given redshift, the metallicity increases with the col- 
umn density, since high Nhi implies deep and massive 
potential wells. In fact, it follows closely the behaviour of 
the rotational velocity as a function of the column density, 
see Fig. 7, since the latter is related to the galaxy luminos- 
ity and metallicity as shown by observations (e.g. Zarit- 
sky et al.1994). The sharp change for Nhi — 10 20 cm -2 
corresponds again to the transition from ionized shells to 
neutral cores. 

Fig. 13 shows directly the redshift evolution of the 
metallicity associated with Nhi = 10 18 cm -2 (dashed 



Fig. 13. The redshift evolution of the average metallicity of 
the halos associated to Nhi = W 18 cm -2 (dashed line) and 
Nhi = 10 22 cm -2 (dot-dashed lines). For damped systems two 
metallicities Z c and Zh are shown (see text). The data points 
are from Pettini et al.(1997) for the zinc metallicity of damped 
Lyman-a systems. 

line) and Nhi = 10 22 cm -2 (dot-dashed lines). For the 
latter case (damped systems) we display both metallic- 
ities Zh and Z c , however the diffuse gas metallicity Zh 
should be the most relevant one. We can see that we ob- 
tain very good agreement with observations by Pettini 
et al.(1997). Moreover, we also obtain a large spread in 
metallicities (in the same way as the velocity distribution 
function was quite extended for a given column density, 
see Fig. 6) which can be seen in Fig. 12: although this pic- 
ture only displays the average Zh, the metallicity disper- 
sion can be estimated from the sharp variation seen near 
Nhi ~ 10 21 cm -2 where one probes different clouds. As 
was suggested by Pettini et al.(1997) in order to explain 
their observations, this is due to the fact that damped sys- 
tems are drawn from a large population of parent galac- 
tic halos which have different star-formation histories and 
physical characteristics. This also provides a good check 
of the validity of our model. 

4-6. Opacity 

Fig. 14 displays the evolution with redshift of the average 
hydrogen opacity thi{z). We can see that our result is 
roughly consistent with observations. In fact, this was al- 
ready implied by Fig. 2 and Fig. 4 which showed that our 
model reproduces the evolution of the column density dis- 
tribution. Indeed, the opacity is directly linked to the lat- 
ter: 

^=( 1+z tef* (44) 
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Fig. 14. The evolution with redshift of the average hydrogen 
opacity thi (lower solid line) and helium opacity T He ii (upper 
solid line and dot-dashed line, see text). The low dashed line 
shows the hydrogen opacity r uni f which would be produced by 
a uniform IGM. The data points are from Press et al.(1993) 
(circles), Zuo & Lu (1993) (filled circles) for hydrogen, and from 
Davidsen et al.(1996) (filled rectangle) and Hogan et al.(1997) 
(cross) for helium. 



where Whi is the equivalent width. As a comparison, the 
dashed line shows the opacity which would be produced 
by a uniform IGM: 



thin and verify Nfjeii /Nhi — 1.8 J ° / J °. Assum- 

J _ ncll > nl _ 912 A 1 228A 

ing that the radiation spectrum is a power-law with in- 
dex -2, Haardt & Madau (1996) get N HeII /N HI ~ 30 
for these clouds. We can sec from Fig. 14 that the helium 
opacity we obtain in this way (upper solid line) is smaller 
than the observations. This discrepancy may be due to a 
change in the slope of the ionizing radiation spectrum. In 
particular, the latter may show strong ionization edges at 
the HI, Hel and Hell ionization frequencies and display a 
step-like profile (see Gnedin & Ostriker 1997, Valageas 
& Silk 1998) which could significantly increase the ra- 
tio N Heii /N hi- Thus, we also show in Fig. 14 the helium 
opacity we get with Nneii /Nhi — 300 (dot-dashed line) . 
We note that, assuming that the Nhi column density dis- 
tribution follows a simple power-law in column density 
and redshift (chosen to be consistent with observations) 
and using N He ii/N H i = 100, Zheng et al.(1998) find that 
at z ~ 3 half of the observed opacity is accounted for by 
clouds with Nhi > 10 12 cm -2 . The helium opacity we 
obtain (dot-dashed line) is higher and consistent with the 
data because i) we use a larger ratio Nneii /Nhi and ii) 
the column density distribution we predict in our model 
extends down to small objects with Nhi < 10 12 cm~ 2 , 
as shown by the lower cutoff Nniiower in (14) and (17). 
Thus, as noticed by Zheng et al.(1998) the large observed 
helium opacity strongly suggests that a significant part of 
the Hell absorption is produced by small density fluctu- 
ations which are below the observational limits for forest 
clouds detected through HI absorption. This population 
is also a natural prediction of our model. 



Tre 2 / m(z) 
m e v a H(z) 



(45) 



Thus, the dumpiness of the distribution of matter de- 
creases the slope of the opacity as a function of redshift, as 
was already noticed by several authors (e.g. Bi & David- 
sen 1997). Indeed, from (24) we expect for a constant UV 
background J21 that 



THI (z)oc(l + z)^n (z) 2 /^-V 
az 



C (A c ) 



(3-2 7 )/(2 7 -l) 



If the clustering is stable, which is a good approximation 
for virialized halos, we obtain: 

THl(z) OC (1 + Z ) 3 ( 6 -T)/( 2 T) = (1 + 2)3-5 for 7 _ L8 ( 46 ) 

which is already very close to observations: Dobrzycki 
& Bechtold (1996) find t H i = 2.6 10~ 3 (1 + z) 3 3 , while 
Tunif a (1 + z) 4 - 5 . We can also obtain the helium II opac- 
ity provided we know the ratio Nneii /Nhi- We used the 
curve Nneii /Nhi calculated by Haardt & Madau (1996) 
for J21 = 0.5. In fact, as was already noticed by Miralda- 
Escude & Ostriker (1992) and Haardt & Madau (1996), 
and is easily checked numerically, most of the opacity 
comes from absorbers with r < 1 which are optically 



4-7. Repartition of matter between different classes of ob- 
jects 

Finally, it is interesting to evaluate the fraction of the mass 
of the universe that is contained in the different popula- 
tions of Lyman-a clouds (it is the same proportion for 
baryonic and non-baryonic matter in our model). This is 
displayed in Fig. 15 as a function of redshift. We can see 
that the mass contained in Lyman limit and damped sys- 
tems (associated with galaxies) increases as time goes on, 
together with structure formation via gravitational clus- 
tering. Moreover, the mass within damped systems is al- 
ways very small, as it corresponds to the small deep cores 
of halos. It increases slightly with redshift because the 
density of all objects is larger at higher redshifts, which 
increases the relative masses of the neutral cores as com- 
pared to the halo masses. Simultaneously, the mass within 
the Lyman-a forest gets larger in the past when most of 
the matter present in the universe is contained in low den- 
sity contrast areas which have not yet virialized (at the 
scale Rd(z)). We must note that even at z = 0, these 
clouds form ~ 30% of the mass of the universe. Thus, at 
any redshift an important part of the mass of the uni- 
verse is contained in small "low density" regions, which 
are not associated with galaxies or luminous matter. This 
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Lyman limit and damped systems occupy a negligible vol- 
ume fraction. 



4-8. Correlation function of various clouds 
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Fig. 15. The evolution with redshift of the fraction of mat- 
ter contained within Lyman-a clouds (for all column densities: 
whole lower dashed area, and for Nhi < 10 11 cm -2 : dashed 
area under the lower solid line), Lyman limit systems (white 
area) and damped systems (upper dashed area). The area 
above the dot-dashed line which decreases sharply for z > 2.5 
is the volume fraction occupied by the Lyman-a forest clouds 

with N HI > 10 11 cm' 2 (Lyman limit and damped systems Fig. 16. The two-point correlation function in velocity space 
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occupy a negligible volume fraction) while the area under this 
curve is the volume fraction occupied by the Lyman-a forest 
clouds with Nhi < 10 11 cm -2 and the very low-density regions 
(or voids) below the threshold Nniio-wer, see (14), which do not 
produce sufficiently large absorption to be observed. 



for clouds with Nhi > 10 cm - , at the redshifts z = (solid 
line), z = 2.5 (dashed line) and z = 5 (dot-dashed line). The 
data points are from Cristiani et al.(1997) (dots) at 1.7 < z < 4 
and from Ulmer (1996) (rectangle) at < z < 1.3. 



matter can only be detected through absorption lines, as 
these small objects may not produce significant features 
in the velocity field. Note that a sizeable proportion of this 
mass is embedded within very low column density regions 
Nhi < 10 11 cm~ 2 . Also, the detected -neutral- hydro- 
gen is only a small (typically 10 -5 ) fraction of the total 
baryon mass for Forest and Lyman-limit objects. Fig. 15 
shows the total baryon mass fraction, including thus the 
dominant, but unobserved, ionised fraction. Since we have 
assumed that the ratio of baryonic to total mass is con- 
stant, the mass fractions in Fig. 15 also correspond to the 
total mass in the Lyman-a clouds. Note that the mass 
within damped Lyman systems is small: they neverthe- 
less stick out prominently in the data because they cor- 
respond to totally neutral hydrogen, a given baryon mass 
fraction being thus (~ 10 5 times) more efficient in produc- 
ing the hydrogen absorption lines. The total mass fraction 
formed by these different populations is close to unity (the 
-negligible- mass fraction which we did not count corre- 
sponds to the luminous galactic cores behind which no 
quasar can possibly be seen). The volume fraction occu- 
pied by Lyman-a forest clouds with Nhi > 10 11 cm -2 
is small for z < 2 since in the highly non-linear universe 
most of the volume consists of very underdense regions. 



Finally, we show in Fig. 16 and Fig. 17 the redshift evo- 
lution of the two-point velocity correlation function . In 
a fashion similar to Cristiani et al.(1997) we relate £ v to 
the spatial correlation function £ Q of Lyman-a clouds by: 

Ur)P(v\r)dr (47) 

with v > 0, r > 0. The index a refers to the fact that £ Q is 
linked to the matter correlation function £ through a bias 
factor b 2 . We use for the conditional probability P(v\r) 
a gaussian of width o~ v centered on Hr = v, where H (z) 
is the Hubble constant and o~ v a characteristic velocity 
dispersion of the considered pair of clouds. The low cutoff 
l m in is simply the sum of both cloud radii l m in = Ri + 
J?2- For small clouds a v -C v so that P(v\r) is narrowly 
peaked at r — v/H, while for large clouds a v > v so that 
the integral is dominated by its lower cutoff, due to the 
divergence at small r of £(r). Thus, we write: 

€v(Nhii,N H i2,v) - Za(N H n,N H i2,r) (48) 

with 

r = M&x(Ri+R2,'^ 7 ^ ] (49) 
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for the velocity correlation function of two lines of column 
densities Njjii and Nhii- We estimate the pair velocity 
dispersion as cr v — (v c i + v c2 )/2, where v C i is the mean 
circular velocity associated with a given Nhh, which was 
described previously in Fig. 7 and Fig. 8. Next, we need 
to obtain the bias parameter b 2 (NHii, Nhi 2 ). Within the 
framework of the scale-invariance of the many-body mat- 
ter correlation functions £ p (i"i, r p ) which led to the 
mass function (4) one can show that at large distances 
r 3> R\ , R 2 j m the highly non- linear regime, the bias char- 
acteristic of two objects factorizes and is a function of the 
sole parameter x introduced in Sect. 2, see Bernardeau & 
Schaeffer (1992). Thus we write: 



Za(N HI1 ,N HI2 ,r) = b(x 1 )b(x 2 )^{r) 



with 



(50) 



x < 1 : b(x) oc x {1 and x > 1 : b(x) oc x (51) 

In our case, a fixed column density Nhi may arise from 
many different clouds, so as for the velocity dispersion we 
shall simply use a mean Xi for each Nhh- Note that the 
previous relation (50) breaks down for r ~ {R\ + R 2 ), but 
in view of the approximations involved (through the use 
of various averages) we shall use it down tor = (Ri + R 2 ) 
where it should still provide a reasonable estimate. Finally, 
since observations usually refer to the mean velocity corre- 
lation function above a given threshold in column density, 
we are led to define: 



&(> N HI ,v) 



iNttl 

dNi dN 2 
n hi *h N 2 



dN 



rj(N) 



(52) 



where T)(Nhi) = d 2 n/ 'dlnNnidz. 

We compare in Fig. 16 the correlation function £ v (> 
Nhi,v) we obtain in this way to observations, for a col- 
umn density threshold Njji > 10 14 cm" 2 as a function 
of velocity. We can see that our predictions agree reason- 
ably well with observations, both at low redshift (data 
from Ulmer 1996) and at high redshift (data from Cris- 
tiani et al.1997). Fig. 17 shows the dependence on the col- 
umn density threshold Nhi of the correlation function 
£„(> Nhi,v) for the redshifts z — 0,2.5 and 5. We re- 
produce the observed increase with Nhi of £ v (> Nhi,v) 
over the range spanned by the data. For small column 
densities, Lyman forest clouds described by the regime 
(1), the correlation function measured for a fixed velocity 
separation increases with Nhi because a higher column 
density corresponds to a deeper potential well, that is a 
higher parameter x, hence a larger bias b(Nni)- Since for 
these low-density clouds the circular velocity is small, as 
seen in Fig. 7 and Fig. 8 (because their virial temperature 
T is lower than T by definition) the spatial separation 




Fig. 17. The dependence on the column density threshold of 
the two-point correlation function £„(> Nhi, 100 km/s) for 
the redshifts z = (solid line), z = 2.5 (dashed line) and 
z — 5 (dot-dashed line). The data points are from Cristiani 
et al.(1997) (disks) for 1.7 < z < 4, Songaila & Cowie (1996) 
(circle) at 2.66 < z < 3.62 and Petitjean & Bergeron (1994) 
(rectangle) at z ~ 2.5. These two latter points are derived from 
CIV metal systems. 



is simply constant: r — v/H(z). Larger column densities 
(Lyman-limit systems) come from deeper potential wells 
with a higher circular velocity so that the bias parameter 
keeps increasing while the separation r becomes influenced 
by the velocity factor a v , hence decreases slowly down to 
Ri+R 2 . Thus the correlation function grows with Nhi- As 
was the case for the mean radius, mass or velocity associ- 
ated to a given column density, £„(> Nhi, v) reaches a flat 
plateau for sufficiently large Lyman-limit systems, when 
different Nhi are drawn from the same population of par- 
ent halos. Finally, there is also a rising part and a higher 
plateau for damped systems, as in the previous studies for 
the velocity or mass. We can notice that there is not a 
deep gap around Nhi ~ 10 20 cm -2 , contrary to these for- 
mer cases, because the correlation function £„(> Nhi,v) 
involves an integral over column densities above a given 
threshold, see (52), which smooths the curves. In fact, 
as can be inferred from Fig. 2 for instance, the correla- 
tion function for lines above 10 19 cm -2 is dominated by 
the contribution of column densities around Nhi ~ 10 21 



5. Open universe: Jlo = 0.3 , A = 

5.1. Column density distribution 

For a low-density universe, f2o = 0.3, A = 0, we can per- 
form the same analysis. Thus, Fig. 18 shows the column 
density distribution at z = 0, 2.5 and 5, while Fig. 19 com- 
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Fig. 18. The contributions (1) corresponding to the Lyman-a 
forest (top left), (2) corresponding to the Lyman limit systems 
(central part) and (2n) corresponding to the damped Lyman 
systems (bottom right), to the column density distribution, at 
z = (solid line), z = 2.5 (dashed line) and z — 5 (dot-dashed 
line) from bottom to top. 




Fig. 19. The total column density distribution (summed over 
all contributions) at z = 3. The data points are from Hu et 
al.(1995) (points) and from Petitjean et al.(1993) (circles). 



Fig. 20. The evolution with redshift of the total column 
density cumulative distribution for Nhi > 10 14 cm~ 2 (up- 
per curve), Nhi > 1.6 10 17 cm -2 (intermediary curve) and 
Nhi > 2 10 20 cm -2 (lower curve). The data points are from 
Bechtold (1994) (crosses), Bahcall et al.(1996) (circles), Gial- 
longo et al.(1996) (filled circles), Kim et al.(1997) (rectangles) 
and Lu et al.(1996) (triangle)for Nhi > W 14 cm -2 ; from Stor- 
rie-Lombardi et al.(1994) (dashed crosses) for Nhi > 1.6 10 17 
cm -2 ; and from Wolfe et al.(1995) (filled circles) and Stor- 
rie-Lombardi et al.(1995) (dashed crosses) for Nhi > 2 10 20 

-2 

cm 

column density distribution we get, and the correspond- 
ing properties of the Lyman-a clouds, come from the ba- 
sic characteristics of the physical model we built to define 
and recognize these absorbers. Of course, the normaliza- 
tion of the power-spectrum plays some role, see Fig. 5 for 
instance, especially for damped systems, but it is not the 
dominant factor for other column densities. Thus, the val- 
ues of J21 (2) we use (see Tab.l) are similar in both cases 
and consistent with observational estimates. 

5.2. Opacity 

Finally, Fig. 21 presents the redshift evolution of the he- 
lium and hydrogen opacities. Of course we obtain results 
very close to the case of a critical universe since the col- 
umn density distributions were already similar. 



pares our predictions with observations at z = 3. Fig. 20 
presents the evolution with redshift of the total column 
density cumulative distribution for Nhi = 10 14 cm~ 2 , 
Nhi = 1.6 10 17 cm~ 2 and N HI = 2 1 20 cm" 2 . We can 
see that we obtain good agreement with the data, simi- 
larly to the previous case of a critical universe. This is not 
really surprising, since as we noticed earlier, provided the 
model-independent properties of the multiplicity function 
discussed in Sect. 2 are satisfied, the main features of the 



6. Conclusion 

In this article we have developed an analytical model for 
the redshift evolution of Lyman-a absorbers. It is based 
on a description of gravitational structures which can han- 
dle objects defined by various density contrast thresholds, 
both negative and positive, assuming that the many-body 
correlation functions are scale-invariant. This allows us to 
go beyond the usual Press-Schechter approach, which is in 
principle restricted to just-virialized halos. This is a key 
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Fig. 21. The evolution with redshift of the average hydrogen 
opacity thi (lower solid line) and helium opacity THdi (upper 
solid line and dot-dashed line). The low dashed line shows the 
hydrogen opacity r uni f which would be produced by a uniform 
IGM. The data points are from Press et al.(1993) (circles), Zuo 
& Lu (1993) (filled circles) for hydrogen, and from Davidsen et 
al.(1996) (filled rectangle) and Hogan et al.(1997) (cross) for 
helium. 

ingredient for modeling Lyman-a absorption lines since 
they arise from a wide variety of objects: from underdense 
clouds to the high-density neutral cores of deep potential 
wells. Once such a tool to count gravitational dark matter 
structures is properly devised, in a consistent way so as 
to keep track of all the mass in the universe while avoid- 
ing double-counting, one still needs to specify a physical 
model in order to associate with each neutral hydrogen 
column density observed on a line of sight the possible 
objects responsible for this feature. 

Similarly to Rees (1986) we consider that large Nri 
lines come from virialized clouds which we associate with 
galactic halos. The external shells ionized by the back- 
ground UV flux lead to Lyman limit lines, while the neu- 
tral cores (protected from this UV flux by self-shielding) 
correspond with damped Lyman systems. A single cloud 
produces a very extended range of possible column den- 
sities, as a function of the impact parameter of the line 
of sight, due to its steep density profile (p oc y 1 for 
our model of spherical clouds). This has important conse- 
quences since it implies that very different Nhi are drawn 
from almost the same population of parent halos. Thus, 
it leads to characteristic features (plateaus and correspon- 
dence between various regimes) in curves such as the mean 
halo rotational velocity or mass as a function of the col- 
umn density. Note that these objects are not defined, 
as is often done, by the sole constant density threshold 
A c ~ 177 as we also take into account a cooling constraint, 
together with the common virialization criterium. 



Moreover, in addition to this population of clouds iden- 
tified with galactic halos, we include a second class of lower 
density objects. These correspond to weak potential wells 
where baryonic density fluctuations are erased over scales 
of the order of the "damping" length Rd- Thus, contrary to 
the former class of clouds, the impact parameter plays no 
significant role and a given object will produce a specific 
column density. These halos, defined by the length scale 
Rd, cover a large range of densities, from very underdense 
regions to density contrasts A ~ 40, and correspond to 
the Lyman forest lines. As discussed in Appendix, the fea- 
tures giving rise to these Lyman forest lines need not be 
real spherical objects and can be treated simply as absorb- 
ing features along the line-of-sight within our framework. 
Thus, although the mean intersecting length Ru of these 
regions with a given line-of-sight is of order Rd (which is 
why we also call Rd the radius of these "objects"), the ac- 
tual extension of such a filament can be much larger. As 
a consequence, the distance at which two separate lines 
of sight can hit the same cloud can be much larger than 
Rd- This explains why the "radius" Rd we obtain is much 
smaller than some lengths obtained from observations of 
close quasar pairs. 

As we noticed above, this wide variety of physical 
objects makes it indeed necessary to use a description 
of gravitational structures which allows one to consider 
in a unified manner a large range of density contrasts 
and scales. For instance, one cannot study Lyman forest 
clouds, nor high column densities at low z, by looking only 
at just- virialized objects, while smoothing the density field 
over the Jeans length prevents one from modelling Lyman- 
limit and damped systems which correspond to very high 
density contrasts at similar scales. 

We also have compared the redshift evolution of the 
column density distribution and the hydrogen and helium 
opacities predicted by our model to observations, for a 
critical CDM universe as well as for an open universe 
S7o = 0.3. In both cases we get a good agreement with the 
data, while the main new parameter is the UV flux J21. 
Indeed, the slope of the column density distribution and 
the relative characteristics associated to various Njji are 
mainly given by the physical model itself we built to iden- 
tify Lyman-a clouds. We can also note that the amplitude 
of the UV flux we need is consistent with observations and 
usual models of structure formation. 

Next, we used the power given by such an analyti- 
cal approach to study the influence of various parameters 
like the normalization of the power-spectrum er 8 , the bary- 
onic density parameter Qj, and the amplitude of the UV 
flux. Thus, we showed how their influence varies according 
to the considered regime (forest, Lyman limit or damped 
systems) which in principle allows one to remove for in- 
stance the degeneracy associated with fi§/J2i when one 
is restricted to the forest and Lyman limit contributions. 
Note however that fih and are not really free param- 
eters, since they are chosen so as to be consistent with 
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an earlier model of galaxy formation. Then, we looked 
at detailed predictions of our model such as the mean 
halo mass or radius associated with a given hydrogen col- 
umn density, as well as the velocity or radius distribution 
themselves for fixed Nhi- This shows that, due to the role 
of the impact parameter, such distributions are very ex- 
tended so that even at z = 2.5, the characteristic halo ro- 
tational velocities associated for instance with Nhi = 10 17 
crn^ 2 cover a large range from ~ 10 km/s up to 500 km/s. 
This also means that one would need very high-resolution 
numerical simulations to take into account the contribu- 
tions produced by all these clouds. Moreover, we showed 
that even within the framework of the "simple" hierarchi- 
cal structure formation scenario the redshift evolution of 
quantities like the maximum halo rotational velocity as- 
signed to a given Nhi displays a qualitatively non-uniform 
behaviour, due to the role of non-gravitational processes 
which introduce additional length or temperature scales. 
Taking advantage of the fact that our model for Lyman-a 
absorbers is part of a broader unified description of the 
structures formed in the universe, including galaxies, we 
used the model of galaxy formation and evolution devel- 
oped earlier (VS II) to obtain the characteristic metallici- 
ties of Lyman- limit and damped systems. We again obtain 
a good agreement with observations, which confirms the 
validity of our approach (both for galaxies and Lyman-a 
clouds !). 

Finally, we considered the redshift evolution of the 
mass and volume fractions formed by the various popula- 
tions of clouds. We note that, contrary to some numerical 
studies, we managed to obtain reasonable agreement with 
the observational column density distribution function to- 
gether with a UV flux J 21 = 0.7 at z = 2.5 consistent with 
observations using a baryonic density parameter = 0.04 
that is close to nucleosynthesis bounds. As we noticed pre- 
viously this could be explained by the very wide range of 
parent clouds which contribute to a given Nhi, so that it 
is difficult for simulations to keep track of all objects (par- 
ticularly weak potential wells) and not to underestimate 
the column density distribution function. Our formalism 
is also very convenient for studying the correlation func- 
tion of Lyman-a clouds. We again obtain results in good 
agreement with observations, both for the redshift evo- 
lution and for the dependence on column density of the 
amplitude of the correlation function. 

Thus, we conclude that the physical picture on which 
our model is based should provide a good description of 
the processes at work in the real universe, since its pre- 
dictions agree with observations for many different quan- 
tities. Moreover, it allows one to get very detailed results 
and keep track of the influence of various processes, while 
building a unified consistent picture of the universe. Of 
course, in order to obtain a simple analytical model we 
had to make some approximations: for instance we con- 
sider spherical clouds (although we tried to correct this 
in a crude way in the regime where it may make a differ- 



ence: at the transition between ionized shells and neutral 
cores) and we did not include the effects of star formation, 
that is likely to be at the origin of the UV flux. Since the 
latter is a key ingredient at all redshifts, it would be inter- 
esting to see whether our model of galaxy formation can 
produce the UV flux needed to match the observational 
constraints on Lyman-a clouds, this will be the subject of 
a forthcoming article (Valageas & Silk 1998). However, the 
amplitude J21 we used is always consistent with observa- 
tional estimates, so that the results obtained in this article 
are quite robust in this respect. Eventually, the detailed 
predictions of our model will have to be checked against 
more precise future data, in order to narrow the range of 
possible physical and cosmological parameters. 
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APPENDIX 

A. Objects versus density fluctuations 

In Sect. 3.1 we have calculated the number of Lyman— a 
forest absorbers along the line-of-sight as if the latter were 
actual distinct objects of size Rd- We show here that there 
is a very simple condition for such a calculation to be 
equivalent to the more direct statistics of the density fluc- 
tuations along the line-of-sight. Indeed, both approaches 
are intimately related and there is a single condition for 
the latter to reduce to the former one. This also shows 
that we take into account all the matter along each line- 
of-sight (we do not restrict ourselves to the highest density 
peaks neglecting low density regions). 

The estimate (2) of the number of objects of size R 
and mass M is based on the statistics of the counts in 
cells of size R. If the fluctuating density field corresponds 



statistics of the counts-in-cells is nothing but the statistics 
of J S(r)d 3 r within the cell. On the other hand, the statis- 
tics of the density field along the line-of-sight, that we are 
interested in here, is the statistics of the smoothed density 
field over the scale Rd, that is the statistics of the aver- 
age Sd = J 8(r)d?r /Vd where Vd is the volume over which 
we average. The probability for the random variable Sd to 
take a given value S at a point r along the line-of sight 
thus has the same expression as the probability to find 
the overdensity 5 within a cell of size Rd, the latter be- 
ing given by the statistics of the counts-in-cells (Balian & 
Schaeffer 1989). It can be written: 



1 dx 
p(S)dS = = xh(x) — 

where x = ^rrrr- The associated column density is 



N HI = 2 n HI R n = 2n lX 2 £{R d f R 



II 
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where Rn is the length along the line-of-sight over which 
the contrast <5 is maintained before the density drops to 
a negligible value. We must recall here that at the small 
scales we consider here, all the matter is within very dense 
spots of negligible volume surrounded by a nearly empty 
space. This is discussed in great detail by Balian & Scha- 
effer (1989) and more recently, in connection with the 
present work, by Valageas & Schaeffer (1997). This equa- 
tion can be usefully compared to (10). Also, the number 
of occurences per redshift interval dz is: 



dii = c^~ — — p(5)d6 



dzR 
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dt dz 1 , , N dx . ^ . 

c———= xh(x) — (A3) 

dz Rn £(R d ) x 



The column density distribution then reads 



d 2 n 



1 1 dt 1 hf \ 



dlnN HI dz) 1 2R// dz £(R d ) 



(A4) 



an expression that is written for a fixed value of Nhi and 
in which x thus depends on R// through (A2). This equa- 
tion again can be usefully compared to (12). In the power- 
law regime h(x) cx x" -2 which is relevant here, it differs 
from the latter by a factor 



±Rd_ 

3i? 7/ 



2 R d 
3R f/ 



(A5) 



to a density p(r), that is to a deviation 5(r) = 



_ p(r)-p 



the 



Provided Rn is governed by the scale Rd, this factor is of 
order unity and our previous estimate in Sect. 3.1 leads to 
the same results as in the present case. 

Next, we can argue that Rn actually is of order Rd- 
Indeed, since density fluctuations (and the density of typ- 
ical matter condensations) increase at smaller scales (fol- 
lowing the increase of the two-point correlation function 
£(i£)) one expect the smallest scale at which fluctuations 
exist (which we called Rd) to dominate. This can also be 
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seen from a simple model where we assume overdensities 
to be in the form of long straight cylinders of radius Rd 
and length L with L 3> Rd- Then, a given filament can 
produce different column densities as the angle 9 between 
its main axis and the line-of-sight varies (along with the 
length Ri I ~ Rd/ sin 0). Indeed, a straightforward calcu- 
lation along the lines of Sect. 3 gives: 

d 3 n \ 1 dt ( N HI N*"" 1 )/ 2 



dlnN HI dlnR/ /dz J L dz \2niR// 



x 1 



_L_\ Rd_ 

Rll) R„ 



1 - 



Rd N 



R 



II 



- 1/a < A6 > 



This shows that the number of lines, for a fixed col- 
umn density Nhi, is dominated by the small values of 
Rj f. the value of the integral of (A6) over Rj/ to get 
d 2 n/dlnNHidz is determined by the lower cutoff Rd- In 
other words, in order to achieve a given Nhi it is more 
likely to look through a higher density region of size Rd 
than to look through a lower density filament almost ex- 
actly along its main axis. This justifies our calculation 
and the equivalence of (A4) with (12). On the other hand, 
we note that the mean intersecting length Ru obtained 
from (A6), that is assuming a straight cylinder, is of or- 
der L (the integral now diverges at large scales). However, 
although long overdense filaments are quite likely to be 
present, a description in terms of a Rayleigh-Levy ran- 
dom walk (with direction jumps over distances of order 
Rd) rather than a long straight cylinder is certainly more 
accurate. In this case, the average intercept is barely twice 
the thickness Rd and all the moments are finite, of order 
Rd- In simple words, even if the distribution of low density 
contrasts is filamentary, it is sufficiently irregular, with di- 
rection jumps over distances Rd, made of N blobs of size 
Rd, but coherent over a distance L ~ N 1 l 2 Rd Rd, so as 
to allow us to describe the distribution of intercepts along 
the line-of-sight as if these features were distinct objects. 
Thus, the description of the number of Lyman forest lines 
in terms of the fluctuating density field and our descrip- 
tion in terms of "objects" of size Rd are equivalent. 

The main difference of such a description in terms of 
filamentary clouds rather than spherical objects is the dis- 
tance at which two different lines-of-sight can cross the 
same cloud. Indeed, whereas the intercept Rii is of or- 
der Rd, the coherence of the absorption features for two 
neighbouring lines-of-sight will be of order L 3> Rd- 
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